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1.  Introduction 


The  most  fundamental  problem  in  the  entire  theory  of  oscil¬ 
lations  is  to  describe  the  motion  of  a  mass  point,  the  tip  mass 

attached  to  a  spring.  Within  the  classical  theory  of  particle 
mechanics,  the  spring  is  regarded  as  massless,  so  that  it  serves 
only  to  transmit  a  force  to  the  tip  mass.  This  force  typically 
depends  on  the  position  and  velocity  of  the  tip  mass  in  perhaps  a 
nonlinear  way.  In  this  case,  the  motion  is  governed  by  an  auto¬ 
nomous  ordinary  differential  equation.  On  the  other  hand,  if  the 
spring  has  mass,  then  its  motion  as  a  continuum  is  coupled  to  that 


of  the  tip  mass.  If  the  spring  has  a  nonlinear  constitutive  equa¬ 
tion,  then  the  analysis  of  the  resulting  motion,  governed  by  par¬ 
tial  differential  equations,  can  be  formidable  indeed. 

this  paper  we-sfrudy  the  motion  of  both  tip  mass  and  spring 
when  the  mass  density  of  the  spring  is  small  and  when  its  consti¬ 
tutive  equation  describes  nonlinearly  elastic  and  viscoelastic 
materials.  Although  these  constitutive  equations  do  not  account 
for  past  history,  -we  nevertheless^,  prove^ythat  in  the  formal  limit 
as  the  spring's  mass  density  goes  to  zero  the  equation  for  the  tip 
mass  is  an  ordinary  differential  equation  for  elastic  springs,  but 
is  generally  not  so  for  viscoelastic  springs.  For  the  latter,  it 
is  typically  an  ordinary-functional  dif f erent^al^aquation  in  which 
the  force  law  for  the  massless  spring  depends  upon  the^^st  histo¬ 
ry  of  the  motion  of  the  tip  mass.  We  furnish  restrictions^**  the 


material  properties  and  initial  conditions  for  the  motion  to  be 
governed  by  an  ordinary  differential  equation  in  the  usual  form. 
Although  these  restrictions  are  quite  special,  it  is  fortunate 


that  they  are  often  encountered  in  practice. 

A  rigorous  asymptotic  analysis  of  the  full  system  of  equa¬ 
tions  for  nonlineariy  viscoeiast ic  springs  is  carried  out  in 

It  shows  that  the  equations  studied  here  constitute  those  for  the 
leading  term  of  the  regular  expansion.  The  dissipative  mechanism 
associated  with  the  viscoelasticity  plays  a  crucial  role  in  this 
analysis.  A  comparable  justification  of  the  asymptotic  status  of 
elastic  springs  is  not  possible.  The  full  asymptotic  analysis 
relies  on  delicate  estimates  for  quasilinear  parabolic  equations. 

The  mathematical  tools  used  in  the  present  paper  are  far  less 
esoteric.  The  novelty  of  our  results  lies  not  in  the  analysis, 
but  in  the  surprising  structure  of  the  equations  and  in  the  meth¬ 
ods  for  manipulating  them.  Nevertheless,  the  finding  of  paradoxi¬ 
cal  behavior  here  would  be  vacuous,  were  it  not  backed  up  by 
rigorous  analysis. 

This  paper  has  two  goals:  (i)  to  describe  the  surprisingly 
complicated  behavior  of  massless  viscoelastic  springs;  (ii)  tc 
serve  as  the  first  step  in  a  development  (pushed  much  further  in 
[1])  of  effective  techniques  for  the  treatment  of  the  dynamics  of 
structures  undergoing  large  motions  and  deformations. 

Notation.  We  denote  partial  derivatives  by  subscripts  and  ordi¬ 
nary  derivatives  by  primes.  To  avoid  ambiguities  we  carefully 
distinguish  between  a  function  and  its  values.  Thus  the  function 
N,  appearing  in  Section  3,  is  formally  presented  as 

(1-1)  [0, 1  ]  •  (-1 .  x.)  3  (x,y)  •— h *  N{x,y )  s  r. 

Then  N  (x,y)  is  its  partial  derivative  with  respect  to  its  first 
argument  evaluated  at  (x,y)  .  We  denote  by  x»-*N(x,y)  or  N(*,y) 


the  function  of  x  obtained  from  N  by  fixing  y.  In  our  nota¬ 
tion  the  term  N(s.w  (s,t))  is  unambiguously  given  by  the  chair. 

s  s 

rule  as 

(1.2)  Nx(s,ws(s,t) )  +  Ny(s,ws(s. t) )wss(s,t) . 

Were  we  to  follow  the  common  practice  of  replacing  the  x  in 
(1.1)  with  s,  then  the  meaning  of  Ng ( s , wg ( s , t ) )  would  no  longer 


cor.f iguraticr.  is  confined  to  a  line.  Since  the  treatment  of  gra¬ 
vity  or  of  other  natural  body  forces  offers  neither  challenge  nor 
insight,  we  take  these  forces  to  be  zero.  Me  may  accordingly  re¬ 
gard  the  spring  as  confined  to  a  straight  smooth  horizontal  groove 
One  end  of  the  spring  is  fixed  at  a  point  taken  to  be  the  origin 
of  coordinates  of  a  line  along  the  groove.  A  particle  of  mass  m 
is  attached  to  the  other  end  of  the  spring. 

The  natural  length  of  the  spring  is  taken  to  be  1.  We  iden¬ 
tify  a  typical  material  point  of  the  spring  with  its  distance  s 
from  the  fixed  end  in  its  natural  state.  In  a  motion  of  the 


spring-mass  system,  at  time  t  the  material  point  s  occupies 
the  position  at  distance  s  +  w(s,t)  from  the  fixed  end.  w(s,t) 
is  the  displacement  of  s  at  time  t.  We  require  that 


(2.1) 


w(s.t)  >  -1,  V  S£  [0,1],  V  t  i  0 


to  ensure  that  the  local  ratio  of  deformed  to  natural  length  r, ever 
be  reduced  to  zero.  The  requirement  that  the  end  s  *  0  be  fixed 


is  expressed  by 

(2.2) 


w(0,t)  ■  0. 


We  impose  the  initial  condition 


(  2 . 3a , b) 


w ( a , 0 )  ■  u(s),  wt(s,0)  ■  *>( a)  for  s  e  [0,1] 


Here  <•>  and  v  are  prescribed,  continuously  differentiable  func¬ 
tions  on  [0.1]  with  <  >’  ( s )  >  -1  for  all  s  -  [0,1]  for  consis¬ 
tency  with  (2.1)  and  with  '.(0)  *  0  for  consistency  with  (2.2). 
Let  -n(s,t)  be  the  (component  of  the)  contact  force  (along 


the  groove)  exerted  on  the  material  of  [s,l]  by  that  of  [0,s) 
at  time  t.  Let  >:  p  ( s)  denote  the  mass  density  of  the  spring  per 
unit  reference  length  at  s.  t-  is  assumed  to  be  twice  contin¬ 
uously  differentiable  and  positive  on  [0,1].  c  is  a  small  posi¬ 
tive  number.  We  assume  that  there  are  no  externally  applied  forces 
acting  on  the  spring-mass  system  (since  their  treatment  would  be 
routine).  If  w  is  sufficiently  regular,  the  linear  momentum 
balance  yields  the  integral  form  of  the  equations  of  motion: 

fl 

(2.4)  -n(s,t)  «  cp(K  )*tt«  ,t)&Z  +  mwtt(l,t), 

J  s 

from  which  we  immediately  deduce  the  differential  equation  of 
motion: 

(2.5)  cp( s)wtt(s,t)  =  ns(s,t) 
and  the  boundary  condition  at  s  *  1: 

(2.6)  mwtt(l,t)  -  -n(l.t) . 

The  properties  of  the  material  at  s  are  specified  by  giving 
the  dependence  of  n(s,t)  on  the  history  of  w  (s,*)  up  to  time 
t  by  means  of  a  suitable  constitutive  functional.  Most  of  our 
interest  will  be  directed  to  viscoelastic  springs  of  differential 
type  1  whose  constitutive  equations  are  expressed  in  terms  of 
given  constitutive  functions 

(2.7)  [0,1]*  ( -1 ,® ) x  (-cd, cd)  a  (x,y,z)  * — ♦  N(x,y,Z)  e  5 
by 

(2.8)  n(s.t)  ■  N(s,w  (s, t ) ,w  (s, t )  )  . 

s  s  t 

(We  systematically  use  the  arguments  introduced  in  (2.7)  to  iden- 


.1  .<  .1  >.i  •*!  i.l  i.(  >i| 


.i . •  .»  -.I  '.t  ut 


tify  the  various  partial  derivatives  of  N.)  For  contrast  we 
shall  discuss  elastic  springs  defined  by  (2.7),  (2.8)  with 

(2.9)  N  =  0. 

z 

To  ensure  that  an  increase  in  extension  (strain)  w  for 

s 

fixed  rate  of  strain  w  ^  be  accompanied  by  a  corresponding 
increase  in  the  tensile  force,  we  require  that 


(2.10) 


Ny  >  °- 


To  ensure  that  a  total  compression  be  accompanied  by  an  infinite 
compressive  force  and  that  an  infinite  extension  be  accompanied  by 
an  infinite  tensile  force  we  require  that 


(2.11)  N(x,y,z) 


as  y 


{-.} 


for  x  e  [0,1] 

and  for  z  bounded, 


To  ensure  that  all  motions  of  viscoelastic  springs  are  dissipative 


we  require  that 

(2.12) 


N  >  const  >  0 
z 


For  simplicity  we  assume  that  N  is  continuously  differentiable 
and  that  N(x,*,*)  is  infinitely  differentiable. 

Assumption  (2.12)  implies  that  z *— »N(x,y,z)  has  an  inverse 


(2.13) 


Z(x,y,n) . 


It  follows  from  (2.7)  and  (2.12)  that 


(2.14) 


from  (2.12)  that 
(2.15) 


Z(x,y,n) 


t®  as  n 


Z  >  0, 
n 


and  from  (2.10)  and  (2.15)  that 


Now  let  n  =  N(y.,y,z)  or,  equivalently,  z  =  2(x,y,n).  Suppose 
that  y — >  -1  or  \  ,  while  n  remains  bounded.  Then  (2.11)  and 

(2.12)  imply  that  z — ►  «  or  -®,  respectively.  Thus  (2.11)  and 
(2.12)  also  imply  that 

(2.17)  Z(x,y,n) - ♦  as  y  - *  {**} 

for  n  bounded . 

The  initial -boundary  value  problem  for  w  is  (2.2),  (2.3), 
(2.5),  (2.6),  (2.8) .  The  reduced  problem  corresponding  to  it, 
obtained  by  setting  c  =  0  in  these  equations,  is 

(2.18)  N(s,ws(s, t ) .wgt(s,t) )  =  N(l,ws(l,t) »w8t(l,t) ) , 

(2.19)  w(O.t)  *  0, 

(2.20)  mwtt(l,t)  +  N(l,ws(l,t),  wst(l,t))  *  0, 

(2.21a,b)  W(s,0)  =  «(s),  Wt(l,0)  a  V’(l). 

In  (2.21)  we  retain  only  those  initial  data  needed  to  ensure  that 
(2.18)  -  (2.21)  has  a  unique  solution  for  small  time.  For  compa¬ 
tibility  of  (2.19)  and  (2.21a,b)  we  require  that  u(0)  *  0.  More¬ 
over,  (2.18)  and  (2.21)  yield 

(2.22)  N(s,«'(s) ,wst(s,0) )  «  N(l,w'(l) ,wgt(l,0) ) , 
which  is  equivalent  to 

(2.23)  wst(s,0)  »  Z(s,«'(s)  ,N(l,w'(l)  ,wst(l,0)  )  )  . 

Since  (2.19)  implies  that 

(2.24)  w  ( 0 , 0 )  *  0, 
we  obtain  from  (2.23)  that 


(2.25) 


wt(s,0)  =  Z(t  )  ,N(l,u'(l)  ,wst(l,0)  )df  . 


There  is  no  reason  to  expect  this  w  (*,0)  to  agree  with  i.  pre¬ 
scribed  in  (2.3b).  In  fact,  setting  s  =  1  in  (2.25)  and  using 


(2.21b)  we  obtain 


(2.26) 


Ml) 


Z(K  .*>'(?)  ,  N(  1 ,«'  ( 1 ) , w  ( l , 0 ) )  )d<  , 


which  we  regard  as  an  equation  for  w  (1,0)  in  terms  of  the  ini- 

S  w 

tial  data  of  (2.21).  (Equation  (2.26)  can  be  uniquely  solved  for 
wgt(l ,0) ,  because  the  right  side  of  (2.26)  is  a  strictly  increas¬ 
ing  function  of  w  (1,0)  that  ranges  over  (-®,a>)  with 
w  (1,0)  .)  This  wgt(l,0)  need  not  agree  with  u'(l)  of  (2.3b). 
Once  wst(l,0)  is  found,  we  can  substitute  it  into  (2.23)  to  pro¬ 
duce  w  (s,0),  which  need  not  agree  with  #>(s)  of  (2.3b). 

&  v 

Remark .  Every  three-dimensional  interpretation  of  our  variables 
shows  that  it  is  reasonable  to  assume  that  N  and  rp  are  each 

proportional  to  the  cross-sectional  area  of  the  spring.  Thus  we 

* 

might  be  led  to  replace  N  with  £N  ,  with  c  interpreted  as 
the  cross-sectional  area.  Such  a  parametrization  would  signifi¬ 
cantly  alter  the  character  of  our  equations.  To  avoid  confusion 
on  this  question,  it  is  helpful  to  regard  the  properties  of  the 
spring,  namely  its  density  and  constitutive  function,  as  fixed  and 
to  parametrize  the  problem  by  the  magnitude  of  the  tip  mass.  (In 
our  analysis  the  small  parameter  K  should  be  interpreted  as  pro¬ 
portional  to  the  ratio  of  the  total  mass  (  [  i.*ds  of  the  spring  to 

J  0 

the  mass  m . ) 


A  simple  physical  argument  reinforces  these  remarks.  Consi¬ 


der  the  equilibrium  of  a  spring  subject  only  to  a  constant  end 
load  n{l,t)  =  '  .  Then  the  eq\\i  librium  version  of  (2.3)  would 

imply  that  n{s,t)  =  *  for  all  s.  Thus  no  introduction  of  a 
small  parameter  c  through  a  constitutive  equation  in  the  form 


(2.27) 


n(s , t )  =  £N  (s,w8(s,t) ,wst(s,t) ) 


could  change  the  value  of  the  tension  n(s,t)  (but  of  course  any 
change  in  X  would  do  so).  If  (2.27)  is  used,  then  the  depen¬ 
dence  of  w  on  f ,  captured  by  an  appropriate  scaling,  would 

S 

have  to  accommdate  the  weakening  of  the  material.  The  scaling 

* 

used  would  have  to  depend  on  N  .  The  adoption  of  (2.27)  in  our 
problem  would  force  us  to  confront  the  dynamical  analog  of  this 
difficulty.  The  analysis  of  [1]  gives  a  precise  standing  to  the 
reduced  problem  (2.18)  -  (2.21);  no  attempt  has  been  made  to  jus¬ 
tify  problems  based  on  (2.27). 


10 


3 .  The  Reduced  Problem  for  Elastic  Springs 

We  begin  our  study  of  the  reduced  problem  in  the  degenerate 
case  of  elastic  springs,  which  have  constitutive  equations  or  tb^ 

form  (2.8)  with  N  independent  of  its  last  argument:  N  *  0. 

The  solutions  of  the  reduced  equations  for  elastic  springs  have 
peculiarities  not  shared  by  those  for  other  materials. 

In  view  of  (2.10),  (2.11),  N(x,‘)  has  an  inverse  (-®,®)  s 

n^-»Y(x,n)  e  (-1,®).  Let  us  assume  that  N(x,0)  =0  so  that 
Y(x,0)  =  0.  Then  the  reduced  equation  (2.14)  subject  to  (2.9)  is 


equivalent  to 
(3.1) 


w  ( s , t )  =  Y ( s , N( 1 , w  ( 1 , t ) ) . 
s  s 


Our  goal  is  to  find  an  equation  for  w(l,*)»  which  governs 
the  motion  of  the  end  mass,  when  w  satisfies  (3.1),  (2.19)  - 

(2.21).  Equation  (2.20),  which  we  rewrite  as 


(3.2) 


mwtt(l,t)  +  N(l,ws(l,t))  =  0, 


is  the  natural  place  to  seek  such  an  equation,  but  is  unsuitable 

as  it  stands  because  w  (l,t)  is  not  related  to  w(l,').  We  now 

s 

obtain  such  a  relation. 

From  (3.1)  and  (2.19)  we  obtain 


(3.3 


w(s , t ) 


YU  ,N(l,w  fl,t)  )  )d*  , 


* 

(3.4)  w  ( 1 ,  t )  «  YU  ,N(l,w  (l,t)  )  )4*  *  J(w  ( 1 ,  t)  )  . 

s  s 

Jo 

The  properties  of  N  show  that  J  increases  strictly  from  -1 
to  a  as  its  argument  increases  from  -1  to  *■.  Thus  (3.4)  is 
equivalent  to 


5X5 


(3.5) 


ws(l,t)  *  J  (w(l,t) ) . 

The  substitution  of  (?.5)  into  (3.2)  yields  an  ordinary  differen¬ 
tial  equation  for  w(l,'): 

(3.6)  mwtt(l,t)  +  N(l»J_1(w(l,t) ) )  =  0. 


Initial  conditions  for  w(l,«)  are  given  by  (2.21a). 

Note  that  (3.4)  implies  that  w(l,t)  >  -1.  If  the  spring  is 
uniform,  i.e.,  if  N  ■  0,  then  J  is  just  the  identity.  The 

A 

properties  of  w(l,*)  are  readily  found  by  a  phase-plane  analysis 
of  (3.6).  The  phase  portrait  consists  of  closed  orbits  around  the 
origin.  By  substituting  (3.5)  into  (3.3)  we  obtain  w  in  terms 
of  w(  1 ,  •  )  : 


(3.7)  w(s,t) 

which  reduces  to 


s 

, N(  1 .  J-1  (w(  1 ,  t )  )  )  )dt  , 

0 


(3.8)  w(s,t)=sw(l,t) 

if  the  spring  is  uniform. 

Note  that  in  general  w  does  not  satisfy  even  (2.21a). 
Equation  (3.7)  shows  that  (2.21a)  holds  if  and  only  if 


(3.9) 


o(s) 


YU  ,N(l,J_1(w(l))))dS. 


0 


From  <3. 7)  we  obtain 


(3.10) 


wt(s,0) 


V>(  1) 


s 

0 


In  the  initial  state  given  by  (3.9),  (3.10)  with  in  ( 1 )  =  0,  the 

spring  is  at  rest  under  the  prescribed  end  displacement  <•»(!).  0 

all  initial  configuration  at  rest,  these,  parameterized  by  <•>{!), 


12 


J"  ✓ 


are  the  only  ones  that  can  be  maintained  with  zero  body  force 
Note  that  for  uniform  springs,  we  find  from  (3.8)  that  (3.9), 

(3.10)  reduce  to 

(3.11)  <o(s)  =  su ( 1 )  ,  wt(s,0)  =  stp(l)  , 
which  have  particularly  simple  Interpretations. 


4 .  The  Reduced  Problem  for  Viscoelastic  Springs 
We  now  assume  that  N  satisfies  (2.12). 

If  the  spring  i^  uniform,  i.e.,  if  N.  =  0,  then  we  i  1  • 

ately  read  off  a  solution  to  (2.18): 

(4.1)  wg ( s , t )  =  wg ( 1 , t ) . 

It  follows  from  (2.19)  that 

(4.2)  w ( s , t )  =  wg(l,t)s,  w( 1 , t )  *  wg(l,t),  wt(l,t)  =  wgt(l,t). 

Substituting  (4.2)  into  (2.20)  we  obtain  the  following  ordinary 
differential  equation  for  the  tip  mass 

(4.3)  mwtt(l,t)  +  N(l,w(l,t),wt(l,t))  =  0. 

In  general,  the  solution  (4.1)  satisfies  the  initial  condi¬ 
tion  (2.21a)  if  and  only  if 

(4.4)  <•>  ( s )  *<•>'(  1 )  s  . 

The  initial  condition  (2.21b)  is  discussed  in  the  paragraph  con¬ 
taining  (2.22).  The  interpretation  of  (4.4)  is  the  same  as  that 
for  (3.9). 

We  now  turn  to  the  construction  of  a  solution  w  of  the 
reduced  problem  when  the  spring  is  not  uniform  or  when  arbitrary 
initial  conditions  of  the  form  (2.21)  are  imposed. 

Let  us  set 

(4.5)  y(t)  «  N(l,ws(l,t) rWst(l,t) ) . 

Then  (2.18)  Implies  that 

(4.6)  wst(s,t)  =  Z(s,ws(s, t) ,i(t) ) . 


We  rewrite  (2.20)  as 


(4.7) 


From  (2.29)  we  get 

(4.8) 


mwtt(l,t)  +  Mt)  *  0. 


w( 1 ,t) 


wg (s , t )ds 

0 


so  that  (4.6)  yields 
(4.9)  wt(l.t) 


Z(s,w#(s,t) ,P(t) )ds. 

0 


We  now  differentiate  (4.9)  with  respect  to  t,  use  (4.7),  rear¬ 
range  the  resulting  equation,  and  use  (4.6)  to  obtain 


(4.10) 


^t(t) 


p(t)/m  +  [  Z  Zdx 
J0  y 


where  the  arguments  of  Z,Zy,Zn  are  s , ( s , t ) ,v ( t ) .  Equations 

(4.6)  and  (4.10)  are  a  pair  of  ordinary  differential  equations  for 

t«-»w  (*,t),i(t).  They  are  subject  to  the  initial  conditions 
s 

(  4 . 1  la  ,  b )  w  (S.0)  =<->{s),  MO)  *  N(1,o'(1)  ,w  (1,0)  ) 

S  Si 


where  w  (1,0)  satisfies  (2.26).  In  a  completely  standard  way, 

S  w 

the  Contraction  Mapping  Principle  implies 

4.12.  Theorem.  There  is  a  number  t+  such  (4,6),  (4.10),  (4.11) 
has  a  unique ,  continuously  differentiable  solution  (0,t+)  ^ 
t~w8(-,t).  Mt)  e  C°(  [0,1] ,  (-1,®)  )xR  . 

Remark.  Condition  (2.12)  ensures  that  (4.9)  can  be  uniquely  solved 
for  v ( t )  : 

(4.13)  i  (t)  =  i  #(wg(- ,t) ,wt(l.t) ) . 

If  we  replace  t  in  (4.6)  and  (4.7)  with  (4.13)  we  obtain  a 
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system  of  two  first-order  ordinary  differential  equations  for 

w_(*,t)  and  w  (l,t)  which  can  be  used  in  place  cf  ( .  C  ■  and 
(4.10)  to  produce  a  theorem  analogous  to  4.12.  Note  that  (  .  7  . 

and  Theorem  4.12  imply  that  w  (1,*)  is  continuously  differen¬ 
tiable  on  (0 , t+) . 

Now  let  us  regard  (4.6)  at  an  ordinary  differential  equation 
for  w  (s,*)  when  v  is  regarded  as  given.  Since  Z(x,#,n)  is 

9 

continuously  differentiable,  the  continuation  theory  for  ordinary 
differential  equations  satisfying  the  Caratheodory  conditions  (cf. 
[3,  Chap.  2 ] , [4 , Sec . 1 . 5 ]  )  imply  that  on  any  interval  [0,t++), 
possibly  containing  [0,t+],  on  which  v  is  bounded  (and  measur¬ 


able),  equation  (4.6)  has  a  unique  absolutely  continuous  solution 
w  (s.*)  satisfying  (4.11a)  as  long  as  w  (s,*)  is  confined  to 

S  5 

(-1,®).  Now  the  parameter  s  appears  in  the  initial  value  prob¬ 
lem  (4.6),  (4.11a)  in  the  initial  datum  u>(s)  and  in  the  depen¬ 
dence  of  Z  on  its  first  argument.  By  the  continuous  dependence 
of  solutions  on  the  data  (cf.  [3,  Theorem  2.4.2])  we  find  that  we 

depends  continuously  on  s  for  t  <  [0,t++)  as  long  as  w  (s,*) 

s 

stays  in  (-1,«).  But  condition  (2.17)  ensures  that  w  (s,*) 

s 

cannot  leave  (-1,®)  on  [0,t++)  as  a  sketch  of  the  slope  field 
( t  ,y)  »-»  Z(s,y,p  (t) )  immediately  shows.  Hence  we  conclude  that 

4.14.  Proposition.  On  any  interval  [0,t*+)  on  which  v  is 

bounded  (and  measurable) .  problem  (4,6),  (4.11a)  has  a  unique 

solution  w_ ( s , * )  with  w  continuous  in  s. 

-  s  -  s  -  — 

It  then  follows  from  (4.10)  that  «■'  must  be  continuous •  on 
such  an  interval.  Let  us  combine  (4.6)  and  (4.7)  to  obtain 


(4.15) 


w>t(s.t)  -  Z(s,w#(s,t),  -  mwtt(l,t)). 


I  uC'ri  'itvrt.  aL’It.  Jl.4*.  I'lilikJ  UAJl  AJ  ililAl  f  *  a  A  A'i  ^iVA4i 


a 


We  can  now  invoke  the  continuous  dependence  of  solutions  on 

the  equation  (cf.  [5,  Section  3.2])  to  obtain 


4.16  Theorem . 


On  any  interval  ■[  0  ,  t  *  )  on  which  i  =  -raw  ( 1 ,  •  ; 

-  XT 


is  bounded  and  measurable,  it  is  continuous .  On  such  an  interval. 


roblem  (4,13),  (4.11a)  has  a  unique  continuously  differentiable 


solution  w  (s,*)  with  the  properties  that  w  (s,t)  e  (-i,a>)  for 
s  s 

0  s  t  <  t++,  w  depends  continuously  on  s,  and  w  (s,t)  depends 
s  s 

continuously  on  the  restriction  of  w  t(l,»)  and  thus  of  w(l,*) 
to  [0,t]  for  o  s  t  <  t++.  A  fortiori ,  w  (l,t)  and  w  (l,t) 

S  —  5  v 

depend  continuously  on  the  past  history  of  w ( 1 , • ) . 

We  now  find  conditions  ensuring  that  t++  **  ®.  Let 


(4.17) 


E(x.y) 


,Y 

N(x,r> , 

-  A 


0)dJ?  . 


Condition  (2.10)  ensures  that  E(x,*)'  is  strictly  convex  and 
(2.11)  that  E(x,y) — as  y — >®.  We  strengthen  (2.11)  by  requir¬ 


ing  that 


(4.18) 


E(x,y) 


®  as  y 


-1 


Let  us  multiply  (2.18)  by  w  (s,t)  and  use  (4.5),  (4,7), 
(4.18)  to  write  the  resulting  equation  as 

Ey(s,ws(s,t) )w8t(s,t) 


(4.19) 


+  |^N(s,ws(s,t)  ,wst(s,t) )  -  N(s,ws(s,t)  ,0)jwst(s,t; 


Mt)wst(s,t)  *  -mwtt(  l  ,t)wgt(s,  t)  . 


Integrating  (4,19)  with  respect  to  s  over  [0,1]  and  using 


The  first  term  In  the  braces  on  the  left  side  of  (4.20)  is  the 
kinetic  energy  of  the  tip  mass  and  the  second  term  is  the  poten¬ 
tial  energy  of  the  spring.  In  the  reduced  problem,  the  spring  has 
no  kinetic  energy.  Since  the  right  side  of  (4.20)  cannot  be  posi¬ 
tive  by  (2.12)  we  obtain 

.1 


(4.21) 


|wt(l,t)2  + 


E  { s ,  w  ( s ,  t ) ) ds  s  C 


Here  and  below,  C  represents  a  positive  constant  independent  of 
t,  depending  only  on  the  data  of  the  problem.  Thus  for  t  e 
(0,t++) , 


(4.22) 
Next 

(4.23) 


|wt(l,t) !  s  C,  |w(l,t) |  s  C(l+t) . 


w(l,t)  >  -1  for  t  in  any  compact 

++ , 


subinterval  of  [0,t  ), 


for  if  not,  there  would  be  a  r  e  (0,t++)  such  that 
(4.24)  w ( 1 , t )  *  1  w8(s,t)dsS“ -1  as  tS r . 

'o 

Thus  w  (*,t)  would  converge  to  -1  in  measure  as  X.S  r .  But 
s 

this  result  contradicts  (4.18)  and  (4.21). 

Now  we  show  that  \i  !  cannot  blow  up  in  finite  time.  If 
Mt) — *®  as  t— *t++,  then  (2.18)  and  (4.5)  would  imply  that  for 


each  s  in  [0,1] 


(4.25)  w  ( s ,  t )  - *  *  or  w  '  as  t  -  ? 

as  a  consequence  of  (2.11),  (2.12).  But  (4.8)  would  then  imply 

that 

(4.26)  w(l,t)  - >  ®  or  wt(l,t)  - ►  ®  as  t  - ►  ® 

in  contradiction  to  (4.22).  If  Mt) — ♦ -®  as  t — ►t++,  then  ana¬ 
logously  for  each  s  in  [0,1] 

(4.27a,b)  wg(s,t)  - ►  -1  or  wst(s,t)  - ►  -a  as  t  - ►  t++. 

Now  (4.27b)  cannot  hold  on  a  set  of  postiive  measure  since  it 
would  imply  by  (4.8)  that  wt(l,t) — >  -®  in  violation  of  (4.22). 

Thus  (4.27a)  would  have  to  hold  for  almost  all  s  in  [0,1] 
implying  by  (4.8)  that  w(l,t) — * -1  as  t — *t++.  This  result  is 
incompatible  with  (4.23)  if  t++  <  ®.  Hence  we  have 

4.29.  Theorem .  Let  (4.18)  hold .  Then  Theorems  4.12  and  4.16 
hold  wl th  t+  =  t++  =  v . 

Theorem  4.16  says  that  w  (l,t)  depends  on  the  past  history 

s 

of  w  (1,*).  If  this  dependence  is  substituted  into  (2.20),  we 
might  well  get  a  functional-dif f erential  equation  for  the  tip  mass. 
Could  the  abstraction  (in  the  form  of  the  contraction  mapping  prin¬ 
ciple)  underlying  the  proof  of  Theorem  4.16  obscure  the  possibili¬ 
ty  that  w#{l,t)  might  depends  only  on  the  present  value  w(l,t)7 
What  are  the  conditions  for  this  to  happen?  I.e.,  what  are  the 
conditions  for  the  motion  of  the  'end  mass  to  be  governed  by  an 
ordinary  differential  equation?  We  address  these  questions  in  the 
rest  of  this  paper.  We  begin  by  studying  an  illuminating  example. 


5 .  Example  Linearly  Vlscot  .astlc  Springs 

We  illustrate  the  results  of  Section  4  with  the  reduced  prob¬ 
lem  for  linearly  viocoelaotic  springs  for  which  :i  ho3  the 


(5.1) 


N(x,y,z)  *  p ( x } z  +  q{x)y 


where  p  and  q  are  positive  and  continuous  on  [0,1].  For  this 
linear  law,  presumably  valid  for  small  motions,  it  is  inappropri¬ 
ate  and  impossible  to  require  that  (2.11)  hold. 


(5.2) 


r(x)  *  q(x) /p(x) . 


Then  we  can  write  (2.18)  as 


5.3)  [er(s)tws(s,t)J  «  p(l)p(s)‘1etr(s,"r(1,]t[er(1,tws(l,t)j 


from  which  we  ultimately  obtain 


>  4 

(5.4a)  w  ( 1 ,  t )  «  f  ( t )  +  aws(l,t)  -  b(s)  er{8)  {,?”t)w^  l  ,n  )dr?ds 

•*0  J0 

where 


(5.4b)  f(t)  -  e  r(s)t[«'(s)  -  p(l)p(s)“V(l )  ]ds, 


(5.4c) 


( 5 . 4d) 


p( 1 )p(s)_1ds. 


b(s)  ■  p(l)p(»)  [r (»)  -  r ( 1 ) ] . 


To  convert  (2.20)  to  an  equation  for  w(l,*)  alone,  we  must 

solve  (5.4)  for  w  (l,t)  in  terms  of  w(l,*).  If  r  is  constant, 

8 

then  the  solution  is  immediate.  Otherwise,  observing  that  the 
integral  over  (0,t)  on  the  right  side  of  (5.4a)  is  a  convolu¬ 
tion,  we  can  accordingly  solve  it  by  the  Laplace  transform  method 


3 


•'.V.'WV.V  ’ 


to  obtain 


(5.5) 


w  (w,t)  =  [  [w(  q  ,  r;  )  -f  (  n  )  ]g(  t-f?  )  d'< 


where  g  is  the  inverse  Laplace  transform  of 

r  f* 

(5.6) 


L  _  b  ( s  )  ds] 

[a  a-r  ( s )  J 


In  order  for  (2.20)  to  be  an  ordinary  differential  equation  for 
w(l,*)#  we  must  be  able  to  express  ws(l,t)  in  terms  of  w(l,t). 
Since  the  solution  (5.5)  is  unique,  this  possibility  can  occur  if 
and  only  if  g  is  the  Dirac  delta,  i.e.,  if  and  only  if  (5.6)  is 
a  constant  function.  This  happens  if  and  only  if  b  =  0,  i.e., 

if  and  only  if  r  is  a  constant.  (The  constancy  of  r  may  be 
regarded  as  a  reflection  of  the  proportionality  of  the  moduli  p 
and  q  to  a  variable  cross-sectional  area  of  the  spring. )  Thus 
if  r  is  not  constant,  then  the  reduced  equation  for  the  tip  mass 
is  that  for  a  massless  viscoelastic  spring  of  memory  type. 

Even  if  r  is  constant,  the  resulting  equation  for  w  (  1  ,  •  '< 
is  not  autonomous  unless  the  transient  f  *  0,  i.e.,  unless  w(l) 

*  aw' ( 1 )  .  We  can  also  determine  conditions  ensuring  that  the  func¬ 
tional  differential  equation  (occurring  when  r  is  not  constant) 
is  autonomous.  A  sufficient  condition  for  f  ■  0  is  that 

(5.7)  h(s)  ■  w'  ( s )  -  p(l)p(s)"V(l)  -  0. 


! 

I 

’  If  r  is  monotone,  (5.7)  is  also  necessary.  For 

t 

|  all  its  derivatives  vanish  at  0: 

1  ,  r  ( 1  ) 


(5.8 


r ( s )  h ( s ) ds  * 


x  h(o (x)  )o'  (x)dx. 


0 


r  { 0 ) 


if  f  ■  0,  then 


k  *  0,1,..., 
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.v-j*. 


jw 


*^nF/W/W TTO*yr7,j^  wy  F,jf  *yfy  F.vw,fj  wwvw  l 
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where  o  is  the  inverse  of  r.  Since  the  monomials  x  are  com¬ 
plete  (in  L2 ,  say),  equation  (5.7)  follows. 

The  development  encompassing  (5.1)  -  (5.4)  can  be  car:  ie-I  _ 

for  nonlinear  constitutive  functions  of  the  form 

(5.9)  N(x,y,z)  =  F(x,y)yz  +  r(x)F(x,y), 

but  such  N's  cannot  satisfy  (2.11).  We  find  that  if  r  is  con¬ 
stant,  then  the  reduced  equation  for  the  end  mass  is  an  ordinary 
differential  equation.  (We  of  course  cannot  apply  the  Laplace 
transform  to  our  nonlinear  equation  to  effect  a  general  solution 
of  the  analog  of  (5.4).) 


'rEF 


6.  Necessary  Conditions  for  the  Reduced  Problem  for  the  End  Mass 
to  be  Governed  by  an  Ordinary  Differential  Equation 
The  results  cf  the  last  two  sections  suggest  two  questions 

about  the  behavior  of  the  reduced  problem  for  nonlinearly  visco¬ 
elastic  springs: 

i)  What  conditions  on  N  ensure  that  the  reduced  problem 
for  the  tip  mass  (generated  by  (2.20))  be  governed  by  an  ordinary 
differential  equation  (possibly  with  a  nonautonomous  transient 
term)  for  all  initial  data  u , u  (satisfying  (2.26)  at  t  =  0)? 

ii)  What  conditions  on  N  and  the  initial  data  ensure  that 
this  problem  is  governed  by  an  autonomous  ordinary  differential 
equation  for  all  values  of  u(l),y(l)?  (These  data  are  those 
appropriate  for  the  motion  of  a  mass  point  on  a  massless  spring.) 

We  obtain  necessary  conditions,  which  can  be  used  to  show  that  such 
ordinary  differential  equations  cannot  occur  for  most  materials. 

We  begin  with  question  (i).  Since  we  have  more  data  to  vary  for 
this  problem,  we  get  much  sharper  restrictions. 

Evaluating  (2.18)  at  t  =  0  and  using  (2.3)  we  obtain 

(6.1)  N(s,w'(s)  ,»'(*))  -  N(l,«/(1)  .v'(l))  , 

which  is  a  compatibility  condition  for  initial  data  (cf.  (2.22)). 

If  (2.20)  is  to  generate  an  ordinary  differential  equation  for 

w(l,  •  ) ,  then  w  (l,t)  and  w  . (l,t)  must  be  functions  of 

s  st 

w(l,t)  and  w  (l,t) .  Indeed,  if  the  order  of  the  equation  is  not 
to  exceed  2,  then  this  dependence  must  have  the  form 

(6.2a)  ws(i,t)  =  j (w( 1 , t ) ,t) , 


(6.2b) 


w.Ji-t)  *  j  (w(  1  ,t)  ,t  )w  ( 1 ,  t)  +  J  (w(l,t),t) 


where 


(6.2c)  (-1.*)  [0.x)  =-  (w ,f)  . - .  J(w,f) 

is  a  function  to  be  determined.  Let  us  set 


(6.3) 


a  *  <•>(!),  /?  *  (P  ( 1 )  . 


By  letting  t — >0  in  (6.2a,b)  we  obtain 

(6.4)  w'(l)  =  j  (a  ,  0 )  ,  v'(l)  =  3w(a,0)'J  +  J?(a,0). 

Now  we  substitute  (6.4)  into  (6.1),  solve  it  for  v'(s),  integrate 
the  resulting  equation  with  respect  to  s  from  0  to  1,  and 
use  (2.19)  to  obtain 


(6.5) 


1 


3  - 


Z(s.o'(s) ,N(1, j{a,0) , Jw(a,0 )B  +  j  (a , 0 ) ) ) ds . 

0 


(Using  (6.5)  we  do  not  have  to  worry  about  how  y>'  is  determined 
from  (6.1).)  Equation  (6.5)  must  hold  identically  for  all  suffi¬ 
ciently  smooth  functions  <•>  satisfying  <>>(0)  *  o  and  <•>'  >  -l 
and  for  all  real  numbers  3.  Thus  the  Gateaux  differential  (first 
variation)  of  the  right  side  of  (6.5)  with  respect  to  «j  in  the 
direction  r?  must  vanish  for  all  r?  with  *7(0)  =  0: 


1 


0 


(ZyJ?'(s)  +  ZntN^jw(a,0) 


+  Nz(Jww(a'0)/J  +  Jw<(a'0))lr){1)}ds- 
Here  and  below  we  understand  the  arguments  of  Z  to  be  those 
shown  in  (6.5)  and  we  use  N*,Ny,...  to  stand  for  the  values  of 
N,N  , ...  at  the  arguments  of  N  shown  in  (6.5).  By  the  funda¬ 
mental  lemma  of  the  Calculus  of  Variations,  we  obtain  from  (6.6) 


(6.7) 


Z  (s.m'(s).N")  =  Z  (l.w'(l)  ,N')  ,  V 


(The  fundaments  1  lemma  also  prod\ices  a  natural  boundary  cond It: r.r. 
from  (6.6).)  The  differentiation  of  (6.5)  with  respect  to  > 


yields 


(6.8) 


1  =  Z  ds  N*1  (a.O).  V  o'.a.fi. 
n  z  w 


Conditions  (2.10)  and  (2.12)  then  imply  that  j  (a,0)  is  every- 

W 

where  positive.  It  follows  that  j  (a,0 )i3  and  therefore  N1 

W 

range  over  (-®,a>)  as  i3  ranges  over  (-&,&).  Thus  (6.7)  implies 
that 


(6.9) 


Zy(s,«'(s)  ,  n)  =  Z  (!.«'( 1)  ,n) ,  V  s,<./,n. 


Since  «'(s)  is  arbitrary  in  (-1,®)  and  since  s  is  arbitrary 
in  [0,1],  we  conclude  from  (6.9)  that 

(6.10)  Zy(x, y,n)  =  Zy ( 1 { 1 )  , n)  ,  V  (x,y,n)  t  -  0 , 1 ) x  ( -1 ,® ) M -® ,® )  • 

(Applying  this  condition  to  (5.1)  we  immediately  obtain  that  q/p 
=  constant,  in  complete  agreement  with  the  findings  of  Section  5.) 
But  condition  (2.17)  implies  that  (6.10)  cannot  hold.  Hence,  our 
answer  to  question  (i)  is  that  there  are  no  materials  satisfying 

(2.10)  -  (2.12)  for  which  the  reduced  problem  for  the  tip  mass  is 
governed  by  an  ordinary  differential  equation  for  all  initial  data . 

We  now  turn  to  question  (ii).  We  assume  that  N(x,0,0)  *  0. 

We  begin  by  studying  (2.18)  -  (2.21)  subject  to  small  initial  data 
of  the  form 


(6.11) 


w(  s  ,  0 )  =  ^(s).  wt(l,0)  »  *7*^(1) 


where  r/  is  a  real  number  small  in  absolute  value.  As  in  Section 
4  we  can  use  standard  results  from  the  theory  of  ordinary  differ- 


ential  equations  to  show  that  (2.18)  -  (2.20),  (6.11)  have  a  solu¬ 

tion  of  the  form 

(6.12)  w(s,t  ;//)  =  n  w!(s,t)  ( //  "  /  2  !  )  w  “  (  s  ,  t )  +  ...  . 

We  get  the  equations  satisfied  by  wr  by  substituting  (6.12)  into 
(2.18)  -  (2.20),  (6.11),  differentiating  these  equations  r  times 

with  respect  to  r? ,  and  then  setting  7?  *  0. 

Let  us  set  Ny(s)  =  Ny(s,0,0),  etc.  Then  w1  satisfies  the 
linearization  of  (2.18),  (6.11): 

(6.13)  N2(s)w*t(s,t)+Ny(s)w*(s,t)  -  Nz( l)w*t( 1, t )+Ny( l)w* ( 1 , t ) , 

(6.14)  w1 (s , 0)  =  (s ) . 

The  analysis  of  Section  5  shows  that  w* ( 1 , t )  has  the  form 
j(w1(l,t))  if  and  only  if 

( 6 . 15a ,  b)  Ny(s)/N2(s)  =r(const.),  ^(1)  =  <^(1)  £nz  ( 1 ) /i*z  (s)J  ds  . 

^  0 


In  this  case  we  readily  find  that 


( 6 . 1 6a , b )  w1(s,t)  =  ?(s)w1(l,t),  =  w1(l,t) 


(6 . 16c) 


r  (s)  s 


Nz(l)/Nz(s) 


J  [R2(l)/Rz(«  )]« 


Now  w  satisfies 


{sz<n 


rt 


§ 

ws«.t)]  +  ert[NyyC:  )\(r-  .t) 


(6.17a) 


+  2Nyz ( K  )w *(*  ,t)wjt(? ,t) 


+  N  (i,  )w 


K  =S 


0, 


i 


(6.17b) 


w  (s , 0)  =  0 . 


Treating  (6.17)  In  a  manner  just  like  that  of  Section  5,  w&  fir;, 
that  w  (l,t)  depends  only  on  wc ( 1 , t )  and  w  (l,t)  (without 
history  terms  or  nonautonomous  terms)  if  and  only  if 
1 


J  {[ 

A 


Nyy(s)r(s)2  -  Nyy(l)?'(l)2jw1(i, 


n 


(6.18)  + 


2^Nyz(s)Ms)2  -  Ny2(l)?'(l)2]w1(lft)wJ(l,t) 

]  wj  { 1 ,  t )  2| 


+  Nzz(s)^(s)2  -  Nzz(l)^(l)2 


ds  =  0 . 


We  let  t — >0  in  (6.18).  Since  0^(1)  and  (1)  are  indepen¬ 
dent,  we  obtain  from  (6.19)  that 

(6.19)  N  ( s  )  T'  ( s  )  2  ,  N  (s)^(s)2,  N  (s)^(s)2  are  constants, 
yy  yz  zz 

Equations  (6.15)  and  (6.19)  are  two  sets  of  necessary  conditions 
for  the  reduced  problem  for  the  tip  mass  to  be  governed  by  an 
ordinary  differential  equation.  It  is  clear  how  to  continue  this 
procedure  ad  infinitum. 

An  alternative  global  approach  to  problem  { i i )  based  upon  the 

formalism  developed  for  problem  (i)  illuminates  some  of  the  issues 

involved,  and  it  can  directly  produce  the  relationship  between 

w  (l,t)  and  w(l,t)  .  In  this  method  we  replace  (6.2a)  with 
s 


(6.20) 


w  (l.t)  «  j(w(l,t)),  J(0)  *  0 


because  we  require  the  resulting  form  of  (2.20)  to  be  autonomous 
The  second  equation  of  (6.20)  is  consistent  with  the  requirement 
that  N ( x , 0 , 0 )  =  0. 


We  adopt  (6.3)  and  allow  <>  to  depend  upon 


(6.21) 


(•>'  =  or  (s,P)  . 


Then  (6.5)  reduces  tc 

1 

i 

(6.22)  .'3  = 


2  (s  ,a.i  (  s  ,a  )  ,N(  1 ,  j  (a  )  ,  j'  (a  )/J  )  )ds  ,  V  a,  3. 


This  equation  restricts  the  functions  Z(*,y,n),.t,  and  j.  If  we 
set  (i  *  0  in  (6.22),  then  we  can  use  (2.10),  (2,11),  (2.15)  to 
solve  the  resulting  equation  for  j  in  terms  of  a  and  t(*,cr). 

We  can  get  alternative  representations  for  j.  Differentiat¬ 
ing  (6.22)  with  respect  to  fi  we  obtain 

.1 


(6.23) 


J'(a)  = 


N 


ZndsJ  ,  V  a,/?. 


Here  the  arguments  of  Nz  and  are  those  of  N  and  Z  in 

(6.22).  We  adhere  to  this  same  convention  below.  Equation  (6.23) 
(or  (6.22))  may  be  regarded  as  an  ordinary  differential  equation 
for  j  parametrized  by  x , Z .  The  right  side  of  (6.23)  is  inde¬ 
pendent  of  3,  so  we  could  set  3=0  there  for  the  purpose  of 
determining  j.  Conditions  (2.12),  (2.15)  imply  that  j'(o)  >  0. 

Thus  J '(a)3  ranges  over  (-od ,<s)  with  3 .  It  follows  that  (6.22] 
.is  equivalent  to 


,1 


(6.24) 

and  to 

(6.25) 


J'(a) 


Z(s,a*(s,a) ,N( 1 , j (a) ,z) )ds,  V  a ,z 


Z(  1 ,  j  (a  )  ,n)  *  J‘  (a ) 


Z(s,at(s,a),n)ds,  V  a,n 


by  the  remarks  leading  to  (2.13).  The  substitution  of  (6.23)  into 
(6.25)  yields 


28 


'  a 


a 


(6.26)  - —  I-lIHHi.?  } - 

=n;  3  ,  D  !  '  »  ,N(1.  j(<<)  ,0) 


Z(s,ar(s,o) ,n)ds 

J  0 


2  { s  ,o  i  (s.«<)  ,N(1-  j(  0)  Ks 

■\ 


This  equation  would  be  an  identity  if  Z =  0,  j(a)  =  a,  r(s,a)  = 
1.  Thus  this  equation  suggests  how  x  must  compensate,  if  pos¬ 
sible,  for  nonuniformity  of  material  response.  We  get  a  collec¬ 
tion  of  necessary  conditions  for  Z(*,y,n)  and  x  to  satisfy  any 
version  of  (6.22)  by  differentiating  it  repeatly  with  respect  to 
a  and  3  (or  equivalent  variables)  and  then  setting  a  =  0  and 
3  equal  to  any  convenient  value.  In  this  process  we  can  generate 
a  power  series  for  j  and  we  recover  results  like  (6.15)  -  (6.19) 


k 

i 

i 

■  7.  Comments 

■ 

;  We  could  easily  treat  problems  in  which  the  parameter 

i 

|  appears  elsewhere,  provided  its  vanishing  does  not  contribute  ^ 

i 

j  further  reduction  of  the  order  of  the  system.  But  our  methods 

t 

J  cannot  handle  the  important  problem  in  which  the  modulus  of 

I 

I  viscosity  N  would  go  to  zero  with  the  parameter  e.  We  could, 

>  z 

\  however,  treat  problems  in  which  the  spring  is  slightly  nonuniform, 

I 

!  with  c  also  serving  as  a  measure  of  nonuniformity.  In  this  case, 

1 

j  the  reduced  problem  would  be  that  for  a  uniform  spring.  (We  would 

;  thereby  lose  at  this  level  the  interesting  phenomena  of  hysteresis 

► 

!  generated  by  nonuniformity.)  Of  course,  to  treat  small  nonuniform- 

[  ity  in  a  truly  satisfactory  manner,  we  should  properly  character- 

» 

j  ize  it  by  an  independent  parameter,  and  then  study  the  asymptotic 

i 

behavior  of  solutions  on  different  regions  of  the  neighborhood  of 
the  origin  in  the  plane  of  these  two  parameters. 

Our  entire  formal  theory  can  be  carried  out  if  (2.8)  is 
replaced  by  a  constitutive  equation  in  which  n(s,t)  depends  upon 
the  history  of  w  (s,*)  up  to  time  t.  The  study  of  the  reduced 
problem,  which  devolves  on  the  analysis  of  (2.20),  can  be  conducted 
as  in  Section  4.  In  particular,  if  the  spring  is  uniform,  then 
(2.20)  for  such  materials  reduces  to 

(7.1)  H(ws(s,- ) )  •  N(w#(l,- ) ) , 

which  admits  a  solution  in  the  form  w  (s,t)  *  w  (l,t).  The 

s  s 

treatment  of  such  solutions  is  identical  to  that  given  in  Sections 
3  and  4.  The  difficulty  with  materials  with  memory  is  that  if  the 
constitutive  equation  lacks  a  term  like  that  of  (2.8)  satisfying 
(2.12),  then  we  cannot  call  upon  the  theory  of  parabolic  equations 
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to  justify  the  expansions  for  the  full  problem  in  [l].  There  is  a 
small  but  illuminating  mathematical  theory  of  cne-aimens i or.al 
problems  for  materials  v.-ith  memory  (cf.  Fenardy,  Hr  us  a  .0.  \ 

for  references)  that  indicates  the  sort  of  difficulties  we  could 
expect  to  encounter.  For  many  such  materials  there  is  a  threshold 
such  that  if  initial  data  are  smaller  than  this  threshold,  then 
the  solution  is  well-behaved  for  all  time,  whereas  if  the  initial 
data  exceed  this  threshold,  then  the  solution  blows  up  in  finite 
time  (i.e.,  the  body  suffers  a  shock).  Thus  we  could  surmise  that 
an  extension  of  the  results  of  Section  4  and  of  [1  ]  to  such  mate¬ 
rials  would  require  at  the  least  some  restrictions  on  the  size  of 
initial  data.  (Hrusa  &  Renardy  [6]  have  studied  a  class  of  mate¬ 
rials  with  memory  that  comes  closer  to  capturing  the  strong  dissi¬ 
pative  mechanism  of  (2.8),  (2.12)  than  those  just  described.) 

Many  rheological  theories  describing  materials  with  memory 
are  based  on  the  modelling  of  behavior  on  the  molecular  level  with 
systems  of  springs  and  dashpots.  It  is  interesting  to  examine  the 
discrete  model  in  order  to  discern  in  microcosm  some  of  the  diffi¬ 
culties  we  have  faced. 

For  simplicity  we  examine  the  linear  system  with  just  two 
degrees  of  freedom  shown  in  Figure  7.2.  k  and  €  denote  spring 
constants,  b  and  c  denote  damping  constants,  and  cm  and  m 
denote  masses.  Let  x  and  w  denote  the  displacements  from  the 
rest  positions  of  the  masses  cm  and  m.  Then  the  equations  of 
motion  of  this  system  are 


(7.3) 


mx  +  (b+c  )y.  -  cw  +  (k+Ox  -  ■'w  =  0, 


(7.4) 


mw  +  c(w-x)  +  *(w-x)  *  0, 


the  superposed  dots  denoting  time  derivatives. 


<7-5>  at 


...  f  c  1  fk+* 

x(t)  exp  -  [EJ^jexp^ 


O  d  ,  * 

c,  -n:[w(x)ex.p  - 


:}*  St1 


1 

T 


which  is  analogous  to  (5.3).  The  solution  of  this  equation  for 
x(t)  in  terms  of  w(t)  is  elementary.  If  k/b  *  l/c,  which 
corresponds  to  the  constancy  of  r,  the  equation  for  w  obtained 
by  substituting  the  representation  for  x  obtained  from(7.5)  into 
(7.4)  is  a  second-order  ordinary  differential  equation.  Otherwise 
it  is  a  second  order  functional-differential  equation,  its  memory 
term  reflecting  the  equivalence  of  this  equation  to  the  third-orde 
system  (7.3),  (7.4)  of  ordinary  differential  equations  for  f.  -  0. 


(In  the  problems  treated  in  Sections  4  -  6 ,  we  reduced  a  partial 
differential  equation,  essentially  equivalent  tc  an  infinite  sys¬ 
tem  cf  ordinary  differential  equations,  to  a  single  seccr:d--o;' ..  :: 

ordinary  (functional-)  differential  equation.)  Note  the  simplifi¬ 
cations  that  follow  from  (7.3),  (7.4)  when  b  *  0  *  c;  these  are 
analogous  to  the  results  of  Section  3.  The  full  asymptotic  analy¬ 
sis  of  (7.3),  (7.4),  or  of  a  corresponding  system  with  many  degrees 
of  freedom,  would  follow  directly  from  the  asymptotic  theory  of 
ordinary  differential  equations. 

These  results  and  those  of  Sections  4-6  raise  the  following 
question:  Suppose  we  are  given  some  initial  conditions  and  solu¬ 

tions  of  an  evolution  equation  in  finite-dimensional  space.  When 
can  we  tell  whether  the  equation  is  an  ordinary  differential  equa¬ 
tion  of  some  fixed  order?  Equivalently,  we  are  given  the  input 
and  output  of  a  black  box.  How  can  we  tell  whether  the  dynamics 
of  the  black  box  is  governed  by  ordinary  differential  equations? 

The  occurrence  of  memory  effects  in  elastic  wave  propagation 
in  nonhomogeneous  media  is  well  known  (cf.  Chen  &  Gurtin  [2],  Reiss 
[7]).  The  nonhomogeneities  reflect  signals  before  they  reach  the 
boundaries.  These  reflections  contribute  to  an  effective  memory 
term,  the  effect  being  pronounced  in  composite  materials  (as  Chen 

&  Gurtin  show).  For  our  problem,  memory  effects  are  not  manifested 
In  elastic  springs,  but  only  in  springs  with  viscous  dissipation. 

In  1911,  Timoshenko  (cf.  [9,  Section  4.9])  developed  a  simple 
technique  to  account  for  the  mass  of  a  linear  spring.  Our  results 
show  how  much  richer  the  theory  becomes  when  viscosity  and  nonlin¬ 
earity  are  allowed  to  intervene. 


Appro v<  .  :  ■  r  ■  .  .  ,-j  ; 

ui^triLut i -:i  ujiii-*ted. 
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